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Abstract. The influence of rotating binary systems on the light 
curves of galactic microlensing events is studied. Three dif- 
ferent rotating binary systems are discussed: a rotating binary 
lens, a rotating binary source, and the earth's motion around 
the sun (parallax effect). The most dramatic effects arise from 
the motion of a binary lens because of the changes of the caus- 
tic structure with time. I discuss when the treatment of a mi- 
crolensing event with a static binary model is appropriate. It 
is shown that additional constraints on the unknown physical 
quantities of the lens system arise from a fit with a rotating bi- 
nary lens as well as from the earth-around-sun motion. For the 
DUO#2 event, a fit with a rotating binary lens is presented. 

Key words: gravitational lensing — dark matter — binaries: 
general — Stars: low-mass, brown dwarfs — planetary systems 
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1. Introduction 

It is a fact that mass objects in a binary system exhibit a rotation 
around their center of mass. However, this motion has mostly 
been neglected in discussions of binary sources and lenses in 
the context of galactic microlensing. Binary motion can play 
a role through a binary lens, a binary source, and the motion 
of the earth (with the observer on it) around the sun (paral- 
lax effect). Griest & Hu ( 1992 ) have presented an examplary 
light curve for a rotating binary source and have claimed that 
such events are rare. As we will see in this paper, however, 
the influence of the rotation effects increases with the event 
timescale. In addition, for the unknown halo population noth- 
ing is known about the distribution of parameters for hypothet- 
ical binary systems. Taking into account the large uncertainties 
in the position of the lens and the relative velocity between 
source, lens, and observer, the rotation effects cannot be ne- 
glected a-priori and fits with static models should be checked 
for consistency. In fact, an event showing the motion of the 
earth around the sun has already been detected (Alcock et. al. 



in the source, the lens, and in the earth-sun system (at the ob- 
server), where the most dramatic effects are caused by a rotat- 
ing binary lens through the motion of the caustics. 

In Sect. 2, a description of the binary motion is given, which 
is needed in the following sections. Section 3 reviews some ba- 
sics of galactic microlensing. Section 4 shows the parametriza- 
tion for rotating binary lens events and some examples for light 
curves. In Sect. 5, estimates are shown which help to decide 
whether the treatment of a binary lens as being static is appro- 
priate. Section 6 discusses rotating binary sources in a simi- 
lar way as for binary lenses. In Sect. 7, I also treat the earth- 
around- sun m otion (parallax effect) which has been n oted b y 
Gould ( |1992| ) and observed by MACHO (Alcock et al. |1995| ). 
In Sect. 8, it is shown how additional information about phys- 
ical parameters follows from the parameters for a rotating bi- 
nary lens and from the parallax effect. For completeness, it is 
also noted that an additional constraint follows from the finite 
source size, if the physical size of the source is known. Sec- 
tion 9 finally presents a fit for the DUO#2 event using a ro- 
tating binary lens, which uses different paramete rs tha n the 
static binary fits already mentioned (Alard et. al 1995 ; Do- 
minik 1997b). The appendix compares the parameters defined 



in Sect. 7 with the treatment of the parallax event by MACHO 
(Alcock et al. |1995|). 



2. Binary motion 

In order to set the notation, I review some properties of the 
dynamics of binary systems (see eg. Landau & Lifshitz 1969 , 
p. 29f.) in this section. Let us consider an object of mass /xi 
at ri and an object of mass fi2 at r2- The Lagrangian of this 
system is given by 



L = ^^Hrl^ + ^^^2r2^ -V{\ri-r2\) 

where V is the gravitational potential 

M1M2 



1995) and the EROS#2 event can be explained by microlens- V (|ri — r2\) — ~G 



ing of an eclipsing binary (an even more special case of a ro- 



ri -r2 



(1) 



(2) 



tating binary source) (Ansari et al. 1995). It is however doubt- 



ful, whethe r the EROS#2 event is due to microlensing at all 
(Paczyriski 1996 ). In this paper I will discuss binary motion 



This is the Kepler problem. 

Let r be the difference vector, R the coordinate of the cen- 
ter of mass, M the total mass and /i the reduced mass given 
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by 




r = 


ri-r2, 


R = 


UlTi + H2r2 



Changing ^ to ^ + 27r corresponds to one revolution, so that the 
period is given by 



M = /ii + ^2 , 

^J■l^J.2 
^J■ = — : — • 

Ml +M2 

With these definitions, the Lagrangian can be written as 

1 , 



(3) 

(4) T = 2n 

(5) 



GM 



(16) 



Since, for n e 2, x{S_ + n27r) — x{^) and y{^ + n27r) = 
(6) yiO one can subtract full periods n from the given time t and 
solve 



'-Mil' 
2 



t' = t-nT = 



-lir 



F(|r|). 



(7) 



GM 



(f-£sin^) 



The Euler-Lagrange equation for R is 

d 



^^dOL 



dt dR, 



i (««') 



forage [0,27r). With 
[a;J = k with k e'^,k < x < k + 1 
(8) and the period T, Eq. ( |l7[ ) can be written as 



so that R = const., i.e. the center of mass moves uniformly. If 
one chooses a coordinate system with the origin at the center 
of mass, one has R = and therefore 

M2 . Ml 



27r 



T 



C-esin^. 



For ^ = at t = tci one has 



ri 



and r2 ~ 



(9) 27r 



t-tp 
T 



T 



Ml + M2 Ml + M2 

As can be seen from the Lagrangian, rit) gives the motion of Po'" ? = Co at ^ = ^i^ C 

to 
(10) 



= C — £ sin C . 
is obtained for 



a particle of mass /i in the gravitational potential 



T 



v(H)^-G^. 



ti - ^(Co -esin^o), 



(17) 



(18) 



(19) 



(20) 



(21) 



For a gravitationally bound system, the trajectory is an el- 
lipse in a plane perpendicular to the angular momentum L, 
where the origin (center of mass) is in a focus of the ellipse. 
Let e be the eccentricity and a the semimajor axis. With polar 
coordinates (r, ip), the trajectory is given by 



so that 

't-ti 



2tt 



T 



t-t 



1 zvr 



(22) 



r(ip) 



1 + £ COS ip 



(11) 

a(l - £) 
--a{l + e) 



+ Co - £ sin^o = ^ - £ sinC 

will yield a ^ G [0, 27r). 

Let r = |r I and v — \r\. The total energy E, which is the 
sum of the kinetic energy T and the potential energy V 

G^M 



where q — a{l — e^). The minimal value is Tmin 
obtained for (p — 0, and the maximal value is /-max 
obtained for ip = tt. 

Therefore, one can parametrize the curve with a parameter j^ 
Cas 



E = T 



V=\,v^- 



is related to the semimajor axis a by 
G^iM 



2a 



riO = a(l-£cosO- 



(12) 



, and with 



The components along the semimajor axis (x-direction) and the 
semiminor axis (y-direction) follow as 



x{^) — a(cosC — £) , 

The time dependence is given by 



t = 



GM 



(C-£sinC): 



(13) 
(14) 



(15) 



From Eq. (23 ) and Eq. 

_ 27r 

f circ — 7^ a 7 

the velocity as a function of r reads 

v{r) = Weil 



2a 
r 



1. 



The maximal velocity is obtained for r^in as 

^inax — '^v ininj — ^cii 



1 



1 



, . „ J 1 • / n\ r^ 1 and the minimal velocity is obtained for Tmax as 

so that t — corresponds to the point (rnun, 0). One sees that -^ 

in general this equation cannot be solved analytically for ^ to [Y 

yield x{t) and y{t) but instead has to be solved numerically, ^mi" = w(fmax) = Vd, 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 
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Pv 



(29) 



Values for pu for different eccentricities e are shown in Table |l]. 
From the virial theorem, one obtains for the expectation values 
of the kinetic and the potential energy the relation 



2 <T>= - <V> 



(30) 



so that one obtains for the radius r and the velocity v the rela- 
tions 



1 1 
r a 



and 



<i;^>= 



GM 47r2 



2^2 



(31) 



(32) 



Table 1. The ratio p^ between the maximal and the minimal 
velocity for different eccentricities e 



e 


Pv 





1.000 


0.01 


1.020 


0.05 


1.105 


0.1 


1.111 


0.2 


1.500 


0.3 


1.857 


0.4 


2.333 


0.5 


3.000 


0.6 


4.000 


0.7 


5.667 


0.8 


9.000 


0.9 


19.000 



3. Some basics of gravitational lensing 

The effect of light bending by a point mass M at the distance 
Dd from the observer and at distance Dds from the source ob- 
ject which is located at a distance Ds from the observer can be 
described by the gravitational lens equation 



y 



\x[ 



(33) 



(e.g. Schneider et al. 1992), where y is a dimensionless coor- 
dinate in the plane perpendicular to the line-of-sight observer- 
lens at the position of the source (source plane) and a; is a di- 
mensionless coordinate in a corresponding plane at the position 
of the lens (lens plane). The physical position of the light ray 
connecting the source and the observer is given by 



r-EX 



(34) 



in the lens plane and 
Ud 



(35) 



In these equations, te denotes the Einstein radius, given by 



r-E 



AGM Dd Dd 



D. 



(36) 



and Tg — ^ fE denotes the projected Einstein radius. For a 
system of N lenses at positions x '^*' with mass fractions nii, 
the lens equation reads 



y 



N 

E 



.{i) 



- X (*) I 



(37) 



\X — X 



Let v± be the velocity of the relative motion between lens, 
source, and observer as measured in the lens plane. If one con- 
siders a coordinate system in which the observer and the source 
are at rest, v± gives the relative motion of the lens. Alterna- 
tively, one can consider a coordinate system in which the ob- 
server and the lens are at rest, so that v± gives the motion of 
the source position as projected onto the lens plane. In either 
case, a characteristic timescale of the motion (and therefore of 
the event) is given by 



iE = ^. 

v± 



(38) 



This definition means that the moving object transverses one 
Einstein radius in the lens plane in the time ^e- 

Let tjjjax denote the time at the closest approach to the line- 
of-sight and Uniin the impact parameter at i^ax in units of the 
Einstein radius te. For a point source and a point-mass lens, 
one obtains for the impact parameter at time t 



<t) - ^<~+WW 



with 



Pit) 



t~U 



tE 



and the amplification is given by (e.g. Paczyhski 1980 

A{u) 



L\Ju^ + 4 



(39) 



(40) 



(41) 



The light curve for an event involving a point source and a 
point-mass lens is therefore described by the 3 parameters tE, 

t-max? anCl liniin- 

4. Rotating binary lenses 

For a rotating binary lens, one needs the projection of the tra- 
jectory onto the lens plane. The orientation of the rotating sys- 
tem relative to the lens plane is given by two angles /3 and 7. 
For /3 = and 7 = 0, a; is chosen along x\,y along xi and the 
angular momentum h is towards the observer (xa-direction). 
The angle /3 describes a rotation of the lens system around x\ 
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and the angle 7 a following rotation of the lens system around 
X2- This means that one has the transformation 



(42) 




with 



7^l 



cos 7 sin P sin 7 cos f3 sin 7 

cos /3 — sin /3 

— sin 7 sin /3 cos 7 cos /3 cos 7 



(43) 



Since z = and the X3 -value is redundant, the transformation 
reduces to 



Xi 
X2 



1 



cos 7 sin /3 sin 7 
cos/3 



(44) 



A rotation around X3 is not considered here, since it can be put 
into the orientation a of the source trajectory. 

Altogether, one needs the following parameters for lensing 
by a rotating binary lens: 

- The point of time if, of the closest approach of the source 
to the center of mass of the lens system, 

- the characteristic time t-E — te/v^, 

- the mimimal projected distance b in the lens plane between 
source and center of mass of the lens system in units of the 
Einstein radius, 

- the angle a between the a; 1 -direction and the direction of 
the projected source trajectory, 

- the mass fraction mi — fii/M, 

- the semimajor axis in units of the Einstein radius p — a/rE, 

- the rotation angle f3, 

- the rotation angle 7, 

- the period T, 

- the eccentricity e, 

- the phase £^0 att = t},. 

Compared with the static binary lens, one needs 5 additional 
parameters. 

The position of a hypothetical object of mass p is therefore, 
using Eqs. (^, ^, and (|4|), 



xiit) 



cos7(cos^(<) — e) + 



sin (3 sin 7 V 1 — £^ sin ^(t) 



X2{t) = pcos/3v 1 — £2 sin^(t) , 

and, with Eq. (ph, the positions of the masses /ii and /i2 are 



(45) 
(46) 



;(i)(t) = (l-mi)a;(t), x'-^\t) = -mix{t) . 



(47) 



From Eq 

t-h 



2tt 



the value of ^ e [0, 2tt) for a given t is given by 

t-U 



T 



T 



tJ^^ 



eo 



esin^o) 
e sin ^0 = C ~ £ sin ^ . 



(48) 



Examples for rotating binary lenses are shown in Figs. |l] 
and g. For both figures, (3 = ^ — e = S,o — 0, and p has been 



chosen as 2x- For Fig. nl the lens model BL for the MACHO 
LMC#1 event has been used (Dominik & Hirshfeld 1996), so 
that tE = 16.27 d, h = 433.18, p = 0.4077, mi = 0.463, 
a = 1.151, h = 0.146, and m,basc = 4.5170. mbasc denotes 
the negative observed magnitude at the unlensed state and / 
the contribution of the source to the total light at unlensed state. 
The rotation period has been chosen as T = 365 d, 100 d, 50 d, 
and 25 d. For Fig. |, the l ens mo del BLO for the OGLE#7 event 
has been used (Dominik [1997^, so that ^e = 80.88 d, h = 
1173.25, p = 1.131, mi = 0.506, a = 2.297, h = 0.048, 
/ — 0.557, and mtaso = —17.5171. The rotation period has 
been chosen as T = 3000 d, 1000 d, 365 d, and 100 d. 

One sees that dramatic effects occur if the period is small, 
especially if additional caustic crossings occur But even for a 
period of 365 days, a deviation from the MACHO LMC#l-fit 
is visible, and for a period of 1000 days a second peak for pa- 
rameters near the OGLE#7-fit occurs. This constellation looks 
a little like the DUO#2 event. A corresponding model is dis- 
cussed in Sect. ^. Note that the rotation period is about 12 times 
larger than the timescale iE for this constellation, nevertheless 
a dramatic effect occurs in the light curve. Note also that the ex- 
ample used by Griest & Hu (1992) for rotating binary sources 
used a period T which is 4 times smaller than ^e- 

5. When is the rotation effect negligible? 

Let us consider a fit for a static binary lens. As discussed by 
Dominik ( 1997a| ), cited as D97a in the following, the rotation 
period T can be estimated using the timescale ^e and distribu- 
tions of the lens position and the velocity v^_. In a similar way, 
one can also obtain probability distributions for the ratio of the 
timescales 

tE 



Rt — 



T 



(49) 



and the ratio of the velocities of the binary motion of the lens 
and the perpendicular motion with respect to the line-of-sight 



i?„ = -^^ . (50) 

Note that irrespective of the eccentricity < v^ >= v^jj.^, and 
the maximal and minimal velocities are of the same order for 
moderate eccentricities as shown in Sect. 2. The projection of 
the velocity of the lens system to a plane perpendicular to the 
line-of-sight may be lower than this, but the velocity is per- 
pendicular to the line-of-sight at least twice a period. Express- 
ing Rt and _R„ in terms of the fit parameters, the lens dis- 
tance X — Dd/Ds, and the dimensionless velocity parameter 
C = v±/vc (where Vc is a characteristic velocity), yields 



Rt — 



and 



tE 



47r y p^ Dg Vc x(l — x) ^ ' 



C / tE 

" 2 Vp^s«ca;(l-x)C 



(51) 



(52) 



This model coincides with that of Alard et. al (1995 ). 
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Fig. 1. Rotating binary lenses with parameters as for the BL-fit of MACHO LMC#1. /3 = 7 = e = ^o = 0. ar = 365 d, b 

T = 100 d, c T = 50 d, d T = 25 d. 



Fig. 2. Rotating binary lenses with parameters as for the BLO-fit of OGLE#7. (3 ^ j ^ e ^ ^o ^ O.aT ^ 3000 d, b T == 1000 d, 
c T = 365 d, d T = 100 d. 
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One sees that 



Rv = 2Trp Rt ■ 



(53) 



Note that there are two problems with these quantities. First, as 
discussed above, the projected trajectory of the binary compo- 
nents are not circles. Second, from fits with static binary mod- 
els, one only gets the projected distance 2% between the ob- 
jects, where p > x in general and p = 2x for circular projected 
trajectories. 

Since Rt, Rv oc ^^ and Rt oc y^TJp^, Ry ex \/l/p, the 
rotation of binary lenses is most important for events with long 
iE and small p^ 

Since Rt and i?„ are of the form 



fc /-i 



G{tE,xX) = Go{tE)[x{l~x)fC 



(54) 



with fc = ? = — i, estimates and probability distributions can 
be derived using the approach presented in D97a. Let H{x) dx 
be the probability for finding x in [x,x + dx] (where H{x) is 
proportional to the mass density p(a;))^, K{() d( be the prob- 
abiUty for finding ( in [(, (^ + dQ, and let T(r, s) be defined 



as 



T{r, s) = / [x{l - x)Y H{x) C K{C) dC dx , 



(55) 



which separates as 

T(r,s)=S(r)W^(s), 
where 

S(r) = f [x{l-x)Y H{x)dx, 



W{s) 



C K{C) dC , 



(57) 



(58) 



if the velocity distribution does not depend on x. Following 
D97a, the expectation values for a quantity G assuming an un- 
known mass distribution are given by 



<G>= Go ^^tn!'^^ = ^0 ^(-1' ^' ■ 



(59) 



Let us adopt the simple galactic halo model of D97a with a 
velocity distribution of 



i^(C) = 2Cexp{-C2} 



(60) 



where Vc has been chosen so that <v^>~ v^, and the mass 
density of halo objects being 



P{r) = Po 



r2 



(61) 



It is implied that one has an effect from binarity. For p ^ 0, 
this effect vanishes. However, if one sees the binarity, the rotation is 
especially important for small p. 

^ H{x) = p{x)/po, where po is an arbitrary characteristic mass 
density, so that at the distance xq, where p{xo) ~ po, one has 
H{xo) = 1. 



where r measures the distance from the Galactic center, Rqc 
is the distance from the sun to the Galactic center, and po is the 
local density at the position of the sun. With the values Ds = 
50 kpc and i?GC = 10 kpc, and the halo being extended up 
to the LMC, located at 82° from the Galactic center as seen 
from the observer, one obtains S(i) = 0.105, S(l) = 0.0407, 
W{-\) = 1.225 and therefore F(-l,-i,-i) = 3.16. Using 
slightly different values for Ds (55 kpc) and i?GC (8.5 kpc), and 
varying the core radius a between and 8 kpc yields estimates 
which differ by about 5 %. 

For the expectation values, one obtains 

1/2 

<Rt> = 1.23 -lO^V"^^^ 



The probability density for 

HR = Rt/<Rt>^ Ry/<Ry> 

is given in D97a as 

4[S(1)]^ 



210km/s 



210km/s 



-1/2 



-1/2 



^ \ 1/2 



(64) 



p{i^r) 



■ exp ■ 



H{x) 

1 S(l) 



KR^x(l-x) 5(i)Ty(-i) 



dx , 



(65) 



(55) and the probability density for Xr = Ig kr is given as 



V'(Ah) = 



[s(i)]= 



4 In 10 



104A« [S(i)M/(- 



H{x) 
)]" J x{l-x) 



■ exp ■ 



2(1) 



lQ^^^x{l-x) S(i)VF( 



dx . 



(66) 



These probability densities are shown in Fig. g, where sym- 
metric intervals around < Rt > or < i?„ > containing proba- 
bilities of 68.3 % and 95.4 % respectively are shown for A^. 
The bounds of these intervals are also shown in Table 0. The 
smallest and the largest value in the 95.4 %-interval differ by a 
factor of about 5. 

For t he bin ary lens fits to MACHO LMC#1 (Dominik & 
Hirshfeld 1996 ), one obtains the values shown in Table BJ Since 
the true semimajor axis a — pr^ is. not yielded by the fit, 
the estimates refer to p = X' which corresponds to a minimal 
value of the period Ty^in, because p > x for any gravitationally 
bound system and T ex p^^^. The timescale t^ ' corresponds 
to the Einstein radius of the smaller mass with mass fraction 
1 



mi. 



tp = is \/l - mi . 



(67) 



One sees that the rotation is not likely to play a dominant 
effect, however a marginal effect may show up. For the wide bi- 
nary models (BA, BAl, BA2), the peak arises from the passage 

(2) 

near the smaller mass on a timescale i^ , so that the influence 
from the binary rotation on the peak will be smaller than esti- 
mated using ts- 
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Table 2. The bounds of symmetric intervals around <Rt> or <i?„> on a logarithmic scale which correspond to probabilities of 
68.3 % and 95.4 % 



AA68.3 


AA95.4 


^0—^-^68.3 


j^Q-^-^es-S 


]^n— ^-^95.4 


]^0^^95-4 


0.1844 


0.3317 


0.654 


1.53 


0.466 


2.15 



Fig. 3. The probability densities a p{kr) (left) and b iP{Xr) 
or <Rv> 



(right) with symmetric 68.3 % and 95.4 % intervals around <Rt> 



Table 3. MACHO LMC#1: Estimates of Rt and Ry for differ- 
ent binary lens models. 





BL 


BLl 


BA 


BAl 


BA2 


tE[d] 


16.27 


17.53 


685 


155 


35.7 


4'' [d] 


— 


— 


17.57 


15.15 


17.72 


<Rt> 


0.054 


0.050 


0.0086 


0.0047 


0.0030 


<Rv> 


0.069 


0.069 


0.13 


0.065 


0.034 



6. Rotating binary sources 

In their discussion of binary sources, Griest & Hu ( 1992| ) have 
also mentioned their rotation. Here I show that the parameters 
for a rotating binary source can be chosen in analogy to the 
rotating binary lens. Let the relative motion of the lens perpen- 
dicular to the source-observer-line projected to the source plane 
be 



y 



(Ions) j'^N ^ / COS a 



smof 






— sma 
cos 5 



6, 



(68) 



where ti, is the point of time of closest approach to the cen- 
ter of mass of the binary source system. The orientation of the 
rotating system relative to the source plane is given by two an- 
gles (3 and 7. For /3 = and 7 = 0, a; is chosen along j/i, y 
along y2 and the angular momentum L is towards the observer 



(i/3-direction). The angle f3 describes a rotation of the source 
system around j/i and the angle 7 a following rotation of the 
source system around j/2- As for the rotating binary lens, one 
has the transformation 



1 



cos 7 sin /3 sin 7 
cos/3 




(69) 



A rotation around y^ need not to be considered here, since it 
can be put into the orientation a of the lens trajectory. 

For lensing of a rotating binary source one needs the fol- 
lowing parameters: 

- The point of time ti, of the closest approach of the lens to 
the center of mass of the source system, 

- the characteristic time ^e ~ 'I'e/v±, 

- the mimimal projected distance b in the source plane be- 
tween lens and center of mass of the source system in units 
of the projected Einstein radius, 

- the angle a between the j/i -direction and the direction of 
the projected lens trajectory, 

- the luminosity offset ratio uj, 

- the mass fraction fhi of source object 1, 

- the semimajor axis in units of the projected Einstein radius 

- the rotation angle /?, 

- the rotation angle 7, 

- the period T, 

- the eccentricity e. 
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- the phase ^q at t = ih. 

Compared with the static binary source, one needs 6 addi- 
tional parameters. 

From Eqs. (^3|), (|l4|), and (|69|), the position of a hypotheti- 
cal object of the reduced mass is given by 



yi{t) = p cos7(cos^(i) -e) + 



y2{t) 



+ sin/3sin7V 1 — £^sin^(i) 
p cos /3 V 1 — e^ sin ^(t) , 



(70) 
(71) 



and the positions of the source objects 1 and 2 are, using 
Eq. (|), 



?/(i)(i) = (l-mi)y(t), y'^^\t)^-miy{t)^ 



(72) 



Note that for ^o = and /3 = 7 = 0, one obtains yi{tb) — 
p(l — e) and 2/2(^6) — 0, so that object 2 is found left from 
object 1 on the j/i-axis, as for the binary lens. The value of 
^ € [0, 27r) for a given t is obtained from (see Eq. (E2b) 



27r 



t-h 
f 






-^0 -esinfo = C-esin^. 



(73) 



The distance of the source objects from the projected lens 
position is given by 



and 



(74) 



(75) 



For a point-mass lens these values can be directly inserted into 
the expression for the magnification A of a point-mass lens (see 
e.g. Paczyriski [1986^ Griest & Hu |1992[) 



A{u) 



i\Jv? + 4 



(76) 



Due to the absence of extended caustics, the effect of a 
rotating binary source and a point-mass lens is less dramatic 
than for a rotating binary lens. Examples are shown in Fig. ^. 
The parameters have been chosen, so that one gets a binary 
source fit for OGLE#5 for T -^ 00 as trans-configuration Q 
This means that ^e = 26.27, a = 1.4512, 4 = 825.719 d, 
h = 0.4624, u = 0.3819, and mbasc = -17.9576. p = 0.8679 
has been chosen so tha^is the distance in the static case and 
m\ = 0.5. /3, 7, e, and ^0 have been chosen as zero and T takes 
the values 100 d, 50 d and 25 d. 



'' The cis-trans-symr netry for binary sources has been discussed by 
Dominik & Hirshfeld (|l996|). 



7. The parallax effect 

7.1. Parameters 

The annual motion of the observer (on the earth) around the sun 
gives another effect of rotating binaries. It has been mentioned 



by Gould (1992) and observed by the MACHO collaboration 



(Alcock et al. 1995). In contrast to the other cases of a rotating 



lens or a rotating source, one knows most of the parameters of 
the binary system: 

- the rotation period T, 

- the semimajor axis a^, 

- the eccentricity e, 

- the point of time when the earth is in perihelion ip. 

One also knows the position of the source of light characterized 

by 

- the longitude (p measured in the ecliptic plane from the per- 
ihelion towards the earth's motion, 

- the latitude x measured from the ecliptic plane towards the 
ecliptic north. 

A fit to an observed light curve for a parallax event involves 
two additional fit parameters: 

- the length of the semimajor axis projected to the lens plane 
measured in Einstein radii p' , 

- a rotation angle t\) in the lens plane describing the relative 
orientation of v^ to the sun-earth system. 

A displacement of the observer's position by (5o is equiv- 
alent to a displacement of the source position projected to the 
lens plane by 



^L = -J^^o = {\-x)5o- 



(77) 



If one chooses x and y in the ecUptic plane and z towards the 
ecliptic north, where the sun is in the origin, positive x is into 
the dkection of the perihelion and positive y is from the perihe- 
hon towards the earth's motion, the motion of the earth is given 
by 



x{m) = ae(cos(C(t))-£), 
y(e(t)) = aeVl-e2sin(e(t)), 
where ^(t) e [0, 27r) can be obtained from 



27r 



T 



T 



C-esin^. 



(78) 
(79) 

(80) 



This motion has to be projected to the lens plane which is to- 
wards the longitude ip and the latitude x as defined before. Let 
the dimensionless coordinates in the lens plane be xi, xi, and 
let xa be a coordinate perpendicular to the lens plane towards 
the observer, so that one gets a right-handed system. If one 
chooses x\ = z,Xi — y and x-^ = —x for (^ = x = 0, 
the angle tp gives a rotation around the xi-axis and the angle x 
a following rotation around the a;2-axis, so that one gets 

1 



Xl 
X2 



rE 



■7^2 



vim) J 



(81) 
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Fig. 4. Rotating binary sources with parameters as for a binary source fit for OGLE#5 (trans-configuration). /3 = j = e = £,0 = 0, 
and mi = 0.5. a T = 100 d, b T = 50 d, c T = 25 d. 



with 



7^2 = 



where 



sm X cos (p 



sm X sm (p 

C0S<y9 



— sm If 
and therefore, with 

P == , 

one obtains 

xi{t) = p' — sinxcos(p(cos^(i) — e) 



— sinxsincpv 1 — e^ sin^(i) 
X2{t) = p — sin93(cos^(t) — e) + 

+ cos(pv 1 — e^ sin^(i) . 
The rotation around X3 by ijj finally gives 

xi{t) — cos V' a?! (t) + sin -0 X2 (0 1 



X2{t) 



- sin -0 xi (t) + cos ip X2 (t) . 



(82) po(i)- 



(83) 



(84) 

(85) 

(86) 
(87) 



Let p be the parameter along the source trajectory and d the 
distance perpendicular to it measured from a line parallel to the 
source trajectory through the origin. By choosing v± along xi, 
I obtain 



p{t) = po{t) + COSIp Xi{t) + Sill'lp X2{t) . 

d{t) = do — sin -0 xi (t) + cos -0 X2 (t) , 



(88) 
(89) 



t-U 



tE 



(90) 



One sees that for t = imax, in general, p{t) ^ and d{t) 7^ dp. 
To avoid a change in the fit parameters when including the par- 
allax (i.e. changing between heliocentric and geocentric coor- 
dinates), it is favourable to use the same fit parameters tmax 
and do (which has been called Umin before) in both cases by 
subtracting the earth-sun distance at imax, which yields for the 
coordinates p{t) towards the direction of the source and d{t) 
perpendicular to it 



Pit) = 
d{t) == 

where 

Po{t) = 



Po(t) + COS-0 (Xi(i) - X2(t,„a; 
+ sin^A {X2{t) - X2(tmax)) , 
do - sinV' (Xi(i) - Xi(imax)) 
+ COS V' (X2(i) - iE2(imax)) , 



:)) 



iE 



and the impact parameter is given by 



u(t) = ^[d{t)Y + WW 



(91) 
(92) 

(93) 
(94) 



The longitude Lp and the latitude x ^6 related to the eclip- 
tic coordinates (3 and A in the following way. The ecliptic co- 
ordinates are geocentric but above a heliocentric system has 
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been used. Therefore, the sun-around-earth motion has to be 
converted to an earth-around-sun motion. It can be seen that 
the vector x of the sun's position measured from the earth is 
transformed into a vector x of the earth's position measured 
from the sun. Since the angular momentum L is an axial vector 
(which means that it does not change its sign under this trans- 
formation), the earth moves around the sun in the same direc- 
tion as the sun moves around the earth in a geocentric system. 
Therefore, one sees that x — P, where the parallax is neglected, 
which does not play a role in determining the position, because 
we deal with distances of the order of 10 kpc. The ecliptical 
length is measured from the vernal equinox along the ecliptic 
in the same sense as the right ascension. Since the sun moves 
towards positive right ascension, A increases with time t. The 
earth's motion around the sun is also in the direction of positive 
A, so that (fi = X + (fie with a constant (/3c, if one neglects the 
earth-sun distance. The sun's position as seen from the earth 
corresponds to A = at vernal equinox, while the earth's posi- 
tion as seen from the sun corresponds to A = tt. If ip^ denotes 
the longitude of the vernal equinox as measured from the peri- 
helion, one obtains ipc = tt + (^^ and therefore ip = X + n + ip-^. 
Inserting the definition of the Einstein radius, Eq. (^6|), into 
the definition of p', Eq. (|3|), yields 



M 



-1/2 



1 



"^RsqDs 



(95) 



a p ex 
= 0.7, 



for X ^ 1. One sees that p' diverges for a: ^ an 
l/^/M. For a; = i, Ds = 8 kpc and M = IMq, p' 
while for M = 10"^ Mq, p' = 23. 

In Fig. pi a light curve where the earth's motion around 
the sun has been considered together with a light curve where 
this motion has been neglected is shown. Since both models 
use the same parameters, the amplification for t = tmax = 
is the same. For this example, parameters which are similar to 
those of the p arallax event found by the MACHO collaboration 
(Alcock et al. 1995 ) have been chosen. The eccentricity of the 
earth's orbit is e = 0.0167, while its rotation period is T = 
365.26 d. With the ecliptical coordinates A = 271° and /? = 
—5°, and ip^ = 77° being the longitude of the vernal equinox 
measured from the perihelion, one obtains ip ~ 2.93 rad and 
X = —0.08 rad. Moreover, I have chosen ijj = 4.14 rad, p' = 



0.2, do = -0.16, tE = 110 d, and t^ 



= to = 0. 



8. Additional constraints on the mass and other physical 
quantities 

For a microlensing event, the Einstein radius te, the lens mass 
M, the lens distance D^ and the transverse velocity v± cannot 
be observed directly in general. Any model for the lens and the 
source involves the timescale ^E; which gives the relation 



TE = tE «_L 



(96) 



so that te can be eliminated and 3 unknown quantities Af , D^ 
and v±, or alternatively p, — M/Mq, x — D^/Ds and v±_ re- 



Fig. 5. Light curves with and without considering the earth's 
motion around the sun, the symmetric curve neglects the earth's 
motion 



main, which are related by the definition of the Einstein radius, 
which reads 



tEWj 



'iGMr; 



© 



Ds V7^*2Kl-aO ■ 



(97) 



Additional constraints may arise from certain models of the 
lens system. In the following, I discuss constraints from the fi- 
nite source size, a rotating binary lens, and the parallax effect. 
Combining two of these allows to determine the lens distance 
(up to a possible degeneracy) and from this value the mass, 
the velocity, and the Einstein radius. Using three or more con- 
straints will overdetermine the problem. However, one should 
note that there are uncertainties in the fit parameters. 

8.1. Using one constraint 

8.1.1. From extended sources 

From the fit of an extended source, one obtains an additional 
constraint if the physical radius of the source r^ is known, 
which may be obtained approximately using the color and the 
absolute magnitude of the source. The parameter _Rsic is the ra- 
tio of the physical radius and the projected Einstein radius r^, 
i.e. 



r„ xr. 



t'e r^ t-E v± 
which is an additional constraint between vj_ and x, i.e. 



v±{x) 



XTs 



tE Rsrc 

Since a; e (0, 1), one obtains a limit for v±: 

v± < 



tERs 



(98) 



(99) 



(100) 



Using the constraint of Eq. (^9|), the lens mass can be written 
as a function of x 



p{x) 



4:GMa)D, 



Rlc 1 



(101) 
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or as a function of v± 

2 
r 

fi{v± 



tEV±K 



4GMqDs i?src(rs - tEV^Rsrc) 



8.1.2. From the observer's motion around the sun 



(102) 



If one takes into account the observer's motion around the sun, 
one has the additional parameter p', which is related to the 
earth-sun distance a® by 



8.2. Using two constraints 

8.2.1. From an extended source and the observer's motion 
around the sun 

Setting equal the expressions for the mass as a function of x for 
the extended source (Eq. (101)) and for the observer's motion 

(113) 



around the sun (Eq. (106)), yields 

^2 



, a®(l-a;) 0^(1 -x) 
P = = —, , 

rE tEV± 

giving a constraint between x and v±, i.e. 

a^il-x) 
v±{x) = -—, . 

tEP 

Since x £ (0, 1), one has 

am 



(103) 



(104) 



(1-X)2 

on 

>Q 



■2 n2 

/2 2 



Since one has 

X 



(114) 



l-x 

for X E (0, 1), only the positive root is an appropriate solution, 

and one obtains 

X a^itsrc 



p'rs 



X . 



l-x 

Solving for x yields the solution 



v±_ < 



tEP' 



1 + X p'n + a©i?. 
(105) and from Eq 



or Eq. (104) one obtains for v\ 



The mass as a function of x follows as (compare Alcock et al. 



v±_ 



affiTs 



tEip'rs +a©i?src) 



(115) 



(116) 



(117) 



1995) 



8.2.2. From an extended source and a rotating binary lens 



x) = 


c2 


2 

a© 
p" 


l-x 

X 




AGMqD, 




d as a 


function ol 


v± 


as 




,,, ^ - 


c2 




tEVi_ 


a© 



4GMqDs p'(ae - tEVj^p') 



8.1.3. From a rotating binary lens 



(106) From the expressions for the mass as a function of x (Eq. ( 101 ) 
and Eq. ( 1 IC )) one obtains the relation 

Ion'' p'^ DsTs 

(107) ^ as a function of x is shown in Fig. pi Note that x e (0, 1). 
The function has zeros for the boundary values x — and 
x = I and a maximum at (I, ^). Therefore, Y is restricted 
to y G (0, ^], which is a consequence from the constraint on 



the Einstein radius, which cannot exceed 



For a rotating binary lens, one obtains from the orbital motion 



r-E,] 



r = 27r 



r, ,3/2 3/2 



GMDs 



(119) 



GM "" ^ y GMqh 

This yields the following pairwise relations: 



p{x) 
vj_{x) 



AnH^p^ 
GMa)T^ 



Q' 



T^c' 



4^6 



1 



10247r4GMQp6 Df x^{l - x)^ 
T^c" I 



(108) For any value of Y , there are two values of x, except for 
y = i^max = -fj- For given x, the mass M, the Einstein ra- 
dius te and the absolute value of the transverse velocity can be 
successively calculated. 

(109) 8 2.3. From the observer's motion around the sun and a 

rotating binary lens 
(110) 



From Eq. ( 1 06 ) and Eq. (110) one obtains 



IGttHeP^ a x{l - x) ■ 
Since x £ (0, 1), v± is restricted to 



(111) x{l-x) 






(120) 



Uj_ > 



AttHeP^D, 



Z as, a. function of x is shown in Fig. R. Note that Z{x) = 
Y(l — x). . The function has zeros for x — Q and x — I 
and a maximum at (i, ^). Since x £ (0, 1), Z is restricted 
(112) to Z G (0, ^]. For any value of Z, there are two values of x, 
except for Z — Z^.^ 



4_ 
27- 
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Fig. 6. r (a;) 



¥ig.7.Z{x) 



8.3. Using three constraints 

With all three constraints, x should be a similar solution to 
X{x), Y{x), and Z{x). Since X{x) yields a unique solution, 
one of the solutions of y(.T) and Z{x) has to be dropped. Note 
that 



X = 



Y 



(121) 



Since the fit parameters and the source radius r^ contain uncer- 
tainties, one may fit for a most-likely value of x. 

9. A fit witli a rotating binary lens for DUO#2 



The DUO#2 event has been reported by Alard et. al (1995) 



where a fit with a strong binary lens is presented. I have inves- 
tigated some more possible fits using a static binary lens and 



a point source (Dominik 1997b). Here 1 show the results of a 



fit with a rotating binary lens with parameters which are com- 
pletely different from the fits with static binaries. 1 have omitted 
one occurence of a data point which appeared twice in the data 
I have received from C. Alard, so that 1 use! 15 data points and 
not 116 data points as in (Alard et al. 1995 1. Moreover, 1 use 



the magnitude values for the fit and not the amplification val- 
ues. The parameters for the fit are shown in Table ^ while the 
resulting light curves are shown in Fig. |[ The x^jn indicates 
that this fit can be marginally accepted. 

Note that the peak after the second caustic crossing is to 
some part due to the rotation and that Rt is 0.07. 

If one adopts a transverse velocity of t^j^ — 30 km/s, one 
obtains, using Eq. (109), a total mass of the lens of M — 
0.025 A/q, i.e. one object with A/i = 5.9 • 10"^ M(^, and one 



Table 4. Fit for DUO#2 using a model with a rotating binary 
lens 



parameter 


value 


te[d] 


6.4 


hid] 


85.40 


P 


1.058 


mi 


0.238 


a [rad] 


2.374 


b 


-0.087 


'^base,blue 


-20.198 


?Tibase,red 


-18.564 


/blue 


0.563 


/red 


0.547 


P [rad] 


0.958 


7 [rad] 


0.358 


r[d] 


92.4 


e 





Co 


-0.046 


Amin 


120.51 


n = # d.o.f. 


100 


\/2xLn-V2n-] 


1.418 


Pix' > Xiin) 


8% 



object with mass AI2 = 0.019 Mq. The Einstein radius be- 
comes TE = 0.11 AU and the semimajor axis a = 0.12 AU. 
However, the lens distance parameter x follows as 0.993 or 
0.007, which seems improbable. Taking 0.993, i.e. the lens 
close to the source, as the more probable configuration, one can 
understand the lens more easily if it is in the bulge population 
rather than in the disk population. However, this value is very 
extreme, a value of 0.9 or 0.95 would have been more plausible, 
but there is some room in the uncertainty of the velocity^ and 
of the fit parameters. Anyway, 1 have tried this model mainly to 
show that there are parameters for a rotating binary model dis- 
tinct from the static ones which produce the light curve, without 
assuming that a physically reasonable model would result. 

As noted in Sect. 5, long timescales iE are favoured for 
showing rotation effects and the timescale for DUO#2 is rather 
short. It is interesting to see how the situation changes if an 
event would have been observed whose light curve shows the 
same shape as that of DUO#2 but whose timescale is 10 times 
larger. For this hypothetical event, the only changes in the 
fit parameters compared to DUO#2 would be an enlargement 
in ^E and T by a factor of 10. For the same transverse ve- 
locity v± = 30 km/s, the Einstei n rad ius r E wo uld be 10 
times larger and, according to Eqs. (109) and (111), the mass 
M and x(l — x) would increase by the same factor, yield- 



^0 



Note that M 
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tory to first order in the eccentricity e. One obtains the earth's 



Fig. 8. DUO #2: Light curve for a rotating binary lens together 
with the data. Light curve for the blue spectral band on the top 
and light curve for the red spectral band on the bottom. 



ing Ml = 0.059 Mq, M2 = 0.19 Mq, te = 1.1 AU, 
a — 1.2 AU, and x = 0.92, i.e. more reasonable values, so 
that the light curve (whose shape is identical to that of DUO#2) 
may well be explained by this rotation effect. This shows that 
it is worth trying fits with rotating binary lenses. 



10. Summary 

Rotating binaries are a reaUty both in the universe and among 
galactic microlensing observations. They are helpful in pro- 
viding additional information about physical parameters of the 
lens system and they may also be used to assign probabilities 
to fits using the knowledge on the distribution of their param- 
eters. The inclusion of the rotation for binary lenses enlarges 
the parameter space and gives room for additional parameter 
degeneracies. It also provides additional shapes of light curves 
through the motion of the caustics. Every fit with a static binary 
should be checked for consistency. 
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Appendix A: Approximation of tlie trajectory to first 
order in e 

Here I show that the ex pressions for parallax light curves given 
by Alcock et. al ( 1995 ) are reproduced by expanding the trajec- 



trajectory in the form (e.g. Montenbruck & Pfleger 1989) 
with 



A{t) = ae( l-ecos(27r^^ 



^(t) = 27r^-:^ + 2esin('27r* ^^ 



With 



A'{t)^Ait) 



T 



1-x 



T 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



the lens plane coordinates xi (t) and X2 (i) follow from Eq. ( pl[ ) 

as 

xi{t) = A'{t)[— sin X cos ip cos £,{t) — 

— sinxsini^sinf (t)] , (A. 5) 

X2{t) = A' {t)[- sin ip cos ^{t) + cos ip sin £,{t)] . (A.6) 

For if = one has 

xi{t) = -A'{t)sinxcos^{t), (A.7) 

X2{t) = A'{t) sinCit), (A.8) 

and therefore 

d{t) = do + A'{t)cos^/jsin^{t) + 

+ A' (t) sin ijj sin X cos ^{t), (A.9) 

p{t) = po{t) — A'{t)cos^psinxcos^{t) + 

+ A'{t)sini:sin^{t). (A. 10) 

For the distance from the origin u{t) one gets 

u^{t) = [d{t)]^ + [p{t)]^ ^ da" + po^ + A'^ sin^ ^ + 
+A'^ sin^ X cos^ ^ + 
+2yl' sin ^ [do cos ip + Po sin ip] + 
+2A' sinxcos^[(io sinip — pq cos?/;] . (A. 11) 

The rotation in the ecliptic plane by an angle ip is equivalent to 
a shift in ^(t). This means that one has 



^(^t) = 2TT^-^+2esin(2TT* *" 



T 



T 



'ip = ^oit)-p>.{A.l2) 



Inserting ^(t) into Eqs. (A.7) and ( |A.8[ ) reveals the expressions 
given in Eqs. (A. 5) and (A.6) where ^(t) has to be replaced by 

The MACHO collaboration (Alcock et al. |1995| ) find that 

u^{t) = ul + uj^{t-tof +a^sin^[n{t~tc)] 

+2asin[n{t — tc)][uj{t — tg) sin9 + uqCosS] + 



+a'^ sin^ P cos'^[n{t - tc)] + 
+2a sin f3 cos[il{t - tc)] ■ 
■[uj{t — to) cos 9 ~ uq sin 9] , 



(A.13) 
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where 

a='^{l-ecos[no{t-t„)]) (A.14) 

V 

and 

n{t - fc) = ^o{t - tc) + 2e sin[f7o(t - tp)] , (A.15) 

where Qq = 27r/r and w = w±/(l — x). 

This expression is equivalent to the derived one with a — 
A', uq = do, ri(t — tc) = £., Po = —oj{t — to), 6 = —tp, 
(fi = rioitc — tp)- Note the sign in po and 0, which is due to the 
fact that I define the lens to be on the right side of the moving 
source, while the MACHO collaboration lets the source move 
in the opposite direction. Further note that their replacement of 
(fi by flo{tc — tp) is an approximation and that tc is the point 
of time where the earth is closest to the sun-source line. Finally 
note that uj — 1/t-E and 

(A.16) 



p'-- 


Lua^ 


Offi 


V 


t^v 
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Figure Ic 
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Figure 2a 
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Figure 2c 
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Figure 3 a 
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Figure 3b 
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Figure 4b 
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Figure 8 
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